Filling-Factor-Dependent Magnetophonon Resonance with Circularly Polarized 
Phonons in Graphene Revealed by High-Field Magneto-Raman Spectroscopy 
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We perform polarization-resolved Raman spectroscopy on graphene in magnetic fields up to 45T. 
This reveals a filling-factor-dependent, multi-component anti-crossing structure of the Raman G 
peak, resulting from magneto-phonon resonances between magneto-excitons and circularly polarized 
F,2 g phonons. This structure is explained with a model of Raman scattering taking into account the 
effects of spatially inhomogeneous carrier densities and strain. This shows that random fluctuations 
of strain-induced pseudo-magnetic fields lead to increased scattering intensity inside the anti-crossing 
gap, consistent with the experiment. 



The magneto-phonon resonances (MPR) effect, ob- 
served in semiconductors when the energy of an opti- 
cal phonon coincides with the inter-Landau level (LL) 
separation [l|, is an important tool for the investiga- 
tion of electron-phonon interactions, especially in two di- 
mensional systems Electron-phonon coupling (EPC) 
in graphene and graphite has been studied for sev- 
eral years 0-0] • The zone-centre, doubly degenerate, 
F,2g phonon strongly interacts with electrons, result- 
ing in renormalization of phonon frequencies and line 
broadening These are tunable by electric and 

magnetic fields, through Fermi-energy shifts and Lan- 
dau quantization. The Raman G peak in graphene is 
predicted to exhibit anti-crossings when the Eo n phonon 
energy matches the separation of two LLs [l44la |. This 
MPR effect can be described as a resonant mixing of 
electronic and lattice excitations into a combined mode, 
leading to a splitting proportional to the EPC[l5j|, and 
was observed in magneto-Raman scattering on SLG on 
the surface of graphite [T^, , non-Bernal stacked mul- 
tilayer graphene on SiC[l9l|. and bulk graphite [20j| . 

Here, we report a polarization-resolved Raman spec- 
troscopy study of MPR in single layer graphene (SLG), 
demonstrating a strong dependence of the MPR line- 
shape on the Raman polarization and carrier density. 
This is explained as a manifestation of MPR between 
electronic magneto-excitons and circularly polarized op- 
tical phonons [15J , combined with the effect of inhomoge- 
neous carrier densities and strain. 

The energy spectrum of SLG in a perpendicular mag- 
netic field B consists of discrete fourfold (spin and val- 
ley) degenerate Landau Levels (LL) with energies E n = 
sign(n) \/l \n\ Hc/Ib , where n = —2, —1, 0, 1, 2, .. is the 



index of LLs in the conduction (n>0) and valence (n<0) 
bands, n = being exactly at the Dirac point (see Refl2ll 
and references therein). The band velocity, c, is the slope 
of the Dirac cone at zero B, and Ib = \JhJeB is the mag- 
netic length. The LL occupancy is characterized by the 
filling factor v = 2iTlgp s , where p s is the carrier density. 
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Fully filled LLs with n = 0,1,... have v = 2,6 
MPR in graphene occurs when the energy of the inter- 
band intcr-LL transitions, — n — > n ± 1, matches that of 
the Fj2 9 phonons, resulting in strongly coupled electron- 
phonon modes, in which the electronic magneto-exciton 
is an anti-symmetric superposition of inter-LL excita- 
tions in each of the two SLG valleys, K and K', while 
this gives symmetric superposition active in far-infrared 
absorption [L5J, |22j. This effect is strongest for the E2 g 
MPR with inter-Landau- level transitions —1—^0 and 
— > 1, at B\ !PR « (25 — 30)T. Since the symmetry 
of such electronic excitation [221] . as well as of the E2 S 
phonons, allows them to be active in Raman scattering, 
MPRs manifest through a fine structure which develops 
in the Raman G peak in the vicinity of the MPR condi- 
Q r = 



tion, 



y/\n\ + i\V2Kc/l B at fl£ 



MPR- 



A specific feature of MPR in SLG, is the v dependence 
of the anticrossing fine structure of coupled electron and 
phonon modes (which can be tuned externally by vary- 
ing the carrier density), as well as a unique possibility 
to resolve the MPR of circularly polarized phonons [15j 
(FigJlJ. The polarization properties of Raman scatter- 
ing involving magnetoexcitons are such that the incom- 
ing and outgoing photons have opposite circular polar- 
izations, with angular momentum transfer ±2S, of which 
±3h is coherently transferred to the lattice, while the re- 
maining angular momentum transfer ^flh goes into elec- 
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FIG. 1. (a) (grey lines) SLG band structure at B=0. In the 
vicinity of MPR between — > 1 electron-hole excitations and 
E2 g phonon, at B l l IPR ~ (25 — 30)T, mode splitting, indicated 
by g, ensues, (b) Calculated [l5t] splitting as a function of filling 
factor, v. Open circles indicate v probed in our experiment 



tronic excitations [22J. In the following, we will refer to 
this angular momentum transfer as a er^-polarized tran- 
sition between the electronic LLs, with opposite signs for 
— n — > n ± 1 inter-LL excitons: cr + for — n — > n + 1 and 
a~ for —n — 1 — > n valcncc-to-conduction band tran- 
sitions. Blocking such transitions by filling (depleting) 
the final states or depleting (filling) the initial states 
with electrons would suppress (promote) coupling be- 
tween magneto-excitons and phonons. Since the a + and 
o~ magneto-excitons are based on different LL pairs, 
changing carrier density in SLG would affect differently 
the size of the anticrossing in the MPR of er T circularly 
polarized phonons, as illustrated in Fig.l. At v = 0, cor- 
responding exactly to a half- filled n = LL, the coupling 
strength of a + - and (T~-polarizcd modes is equal, caus- 
ing the G peak to split equally for a + and o~ phonon 
polarizations For < v < 2, corresponding to a par- 
tially filled n = LL, the — f — > transition becomes 
partly blocked, while the — > +1 transition is promoted, 
giving rise to different splittings in the fine structure of 
<7 + -polarized and o~ polarized modes. For 2 < v < 6, 
n = LL is full, leaving no space for MPR with the a~ 
phonon, whereas n = +1 LL is only partly filled, result- 
ing in a MPR-induced fine structure in the a + phonon 
lincshape (with maximum splitting at v = 2). Thus, the 
EPC at the MPR resolves energies of the er^-polarized 
lattice excitations [la, [22|. Note that, in order to probe 
the direct B influence on the motion of lattice nuclei 
would require unattainably strong B^ 10 4 T or higher, 
as can be estimated from the carbon-atom-to-electron 
mass ratio. Finally, at high v > 6, both — f — > and 
— > 1 transitions are blocked, completely suppressing 
the MPR- induced fine structure for both modes. 



inorganic contaminants from the surface. The foil is then 
loaded in a 4 inch quartz tube, and annealed at 1000°C 
for 30mins with an H2 flow of 20 cubic centimeters per 
minute (seem) and a pressure of 200mTorr. Such an- 
nealing process increases Cu grain size and removes na- 
tive oxide layer. For the growth of graphene films, a gas 
mixture of H2 and CH4 (flow rates of 40 and 20sccm re- 
spectively) is injected in the CVD chamber for 30mins 
while maintaining the pressure at 600mTorr and temper- 
ature at 1000°C. Finally, the sample is rapidly cooled 
to room temperature under a hydrogen flow at a pres- 
sure of 200mTorr. Such procedure results in high qual- 
ity SLG on both sides of the foil, as shown by Raman 
spectroscopy [1, [25|] • A layer of Poly(methyl methacry- 
late) (PMMA) is then spun on one of the surfaces of the 
Cu foil. The graphene layer on the other side of Cu is 
removed with an oxygen plasma at a pressure of 20mTorr 
and a power of 10W for 30s. The PMMA/graphene mem- 
brane is then released by dissolving the Cu in a 0.2M 
aqueous solution of ammonium persulphate. The remain- 
ing PMMA/graphene layer is cleaned several times by DI 
water to remove residual salt, and transferred onto the 
Si02/Si substrates. After the sample is dried in ambi- 
ent conditions, the PMMA is removed in acetone, leav- 



We measure a SLG on Si/Si02 prepared as follows. 
Large area graphene films are grown on copper foils by 
chemical vapor deposition (CVD) 23, 24j]. A^25/^m thick 
piece of Cu foil is initially cleaned with acetone, deion- 
ized (DI) water, and isopropanol to remove organic and 



ing the graphene adhered to Si02/Si[24|. Using the de- 
pendence of graphene Raman spectra upon doping, the 
carrier density is estimated by combining the intensity 
and area ratio of the G and 2D peaks, I(2D)/I(G) and 
A(2D)/A(G), with the position and Full Width at Half 
Maximum of thesepeaks, Pos(G), Pos(2D), FWHM(G), 
FWHM(2D)@,[1, [3. Under ambient laboratory condi- 
tions, the SLG sample is initially p-doped with a car- 
rier concentration^ x 10 12 cm~ 2 . By adjusting anneal- 
ing parameters (typically 8 hours in 90%/ 10% Ar/H2 at- 
mosphere at temperatures up to 220°C) and degassing 
in< 10 _4 mbar vacuum, the sample can be rendered n- 
type, with an electron density^ 2 x 10 12 cm~ 2 . Exposing 
the sample to low-pressure (^lmbar or less) 7V2 atmo- 
sphere reduces the carrier density to^ 0.4 x 10 12 cm~ 2 . 
The exposure to ambient pressure N2 gas or air re- 
sults in hole doping, restoring p-type doping with carrier 
densities^ (5-6) x 10 12 cm~ 2 . If the sample is left in ambi- 
ent air for a long period of time, it continues to experience 
p-doping reaching a carrier concentration^ 10 13 cm -2 in 
several weeks. However, this doping is inhomogeneous, 
with a 10 to 20% variation, as shown by Raman mapping. 

B dependent Raman measurements arc performed in 
a quasi-backscattering configuration at 300K for B up 
to 45T using a high-field magneto-optical insert, as for 
Refl20l. The combination of a linear polarizer and a 
A/4 plate, circularly polarizes the incident and scattered 
light. The excitation beam (532nm) is focused to a 
spot<10/im, with a power^4mW. The circularly polar- 
ized cr^ excitation, probed using the cr ± /(T =F configura- 
tions of out/in polarisations, is achieved by reversing the 
B polarity. In this notation, the first (second) symbol de- 
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FIG. 2. (a,b) Unpolarized magneto-Raman spectra at carrier density corresponding to v > 6. The inset compares 20T a + - 
and a~ circular polarized spectra, (c-f) Circular-polarized magneto-Raman spectra at (c,d) 2 < v < 6 and (e,f) v < 2 for 
B>15T. The dashed lines in (d,f) are guides to the eye for the B dependence of electron-phonon coupled modes at the (0 — > 1) 
MPR anticrossing. The blue curves in (d) are additional spectral components observed in the close vicinity of the intersection 
between a (0 — > 1) magneto-exciton and the E2 9 phonon energy 



fines the polarization of the excitation (scattered) light. 
The polarization efficiency, i.e. the ratio of light trans- 
mitted through two circular polarizers with parallel or 
crossed helicities, is~90% for the incident light, while it 
is~75% for the scattered light in the spectral range of 
interest^l450-1700cm _1 . The resolution is~1.9cm _1 . 

Figs l2l3l show the polarization- and B dependence of 
the G peak at different carrier densities. The Raman 
spectra can be classified in three categories, as a function 
of v. At high carrier density, corresponding to \v\ > 6, 
the G peak does not reveal MPR-induced splitting nor 
any polarization dependence, Figsl2^,,b. As the filling 
factor decreases (y < 6), a significant change in the spec- 
tra is observed. FigsJ2t,d plot spectra measured on the 
sample annealed, degassed and kept in< I0 _4 mbar vac- 
uum. The G peak exhibits a strong, anticrossing-like, 
splitting reaching^ 150cm _1 (~ 20meV) for B=25T. 
The electron-phonon coupled modes appear only in the 
a and a + /<7~ geometry, while the G peak does not 
split nor shift in a~ and er~ / a + polarizations. This in- 
dicates n-type doping with estimated carrier density^ 
2.0 x I0 12 cm~ 2 . The observed behavior reveals the MPR 
at 2 < v < 6. Indeed, the condition v = 6 is met at 
B~10T, with a coupling of — > 1 magnctoexcitons and 
E,2g phonons at higher B. As the carrier density further 



decreases, so that v < 2, the G peak splitting changes 
(FigsJ2k,f). In contrast to the 2 < v < 6 case, the cou- 
pled modes now appear in both a + and a~ polarizations. 
Moreover, the spectra measured at B > 30T reveal that 
the coupled-mode consists of two Lorenzians. Their rel- 
ative intensities depend on the helicity of the incoming 
light. We assign these two components to the polar- 
ized modes originating from the coupling of Fi2g phonons 
with — > 1 and —1—^0 magncto-excitons. We assign 
the fact that polarized modes do not completely dis- 
appear in the measurements with opposite helicity to the 
imperfect polarization selectivity of our setup, combined 
with sample inhomogeneity, as discussed below. 

For the strongest anticrossing, at B^p, the observed 
behavior can be described using the resonance approxi- 
mation formula [l5| for pairs of strongly coupled fi± cir- 
cularly polarized modes: 

.g ff+ (2<^<6)= 5o -^V6-^, ff(T -(2 <i/<6) =0, 

IB 

a . /3 / Ar 

g a ±{0 <u <2) = g —V2±v, go = y -y ^"7o- 
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FIG. 3. Peak position, FWHM and normalized intensity of 
the G-peak Lorenzian components as a function of magnetic 
field. Open black circles represent the central E2 9 phonon 
line component. Filled red and grey circles are the electron- 
phonon coupled modes. Dashed lines plot EqsQ] Blue squares 
show additional components observed in close vicinity of the 
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Here, fl± is the frequency of the upper (+) or lower (-) 
coupled mode, ilr ~ 1581cm -1 is the bare F,2 g phonon 
energy, Q = V^^c/Ib, a (= 2.46A) is the graphene lat- 
tice constant, go describes the EPC at zero B in terms 
of the nearest neighbor hopping integral 70, and the di- 
mensionless coupling constant Ar 0, El ■ 

The set of equations ([T]) describe the energies of the two 
branches of mixed phonon and magneto-exciton modes 
for an idealized non-strained, homogeneously doped SLG, 
also neglecting broadening of the bare excitations due to 
inelastic and elastic scattering. The detailed description 
of MPR and its manifestation in the Raman spectrum of 
realistic samples should include broadening of the opti- 
cal phonon, 7r, and of the magneto-exciton, j e , as well 
as possible inhomogeneity of doping and fluctuations of 
local strain. In a sample with a spatially varying car- 
rier density, g a and g" become coordinate-dependent, 
which would lead to inhomogencous broadening of the 
MPR fine structure. Inhomogeneous strain in SLG gen- 
erates a pseudo-magnetic field. B str , with opposite sign 
in the K and K' valleys 0, [23, El, giving a total field 
B K = B + B str at K and B K , = B - B str at K'. In 
strained graphene nanobubbles, it was reported that B str 
can reach 200-400T[29|. In nominally unstrained SLG, 
the residual strain caused by nm-scale height variations 



and defects can induce random B str >10T[30|. Locally, 
this splits a simple anticrossing into mixing and splitting 
of each <T ± -polarized phonon and two magneto-excitons, 
now distinguished by different LL energies in the two 
valleys, leading to an even more complicated fine struc- 
ture. Also, the quantum efficiency (i.e. the probability 
of a scattering event per incoming photon) of — > 1 
and — 1 — > magneto-excitons in Raman scattering}^ 
is much lower than the quantum efficiency of the T-point 
optical phonon, manifested in the G peak area, A(G)[3l|. 

All these factors can be taken into account in the for- 
mula for the spectral weight of Raman scatterin g at MPR 
derived using the random phase approximation [32|: 



A(G) 
2tt 
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where angle brackets stand for the averaging over the 
spatial fluctuations of a strain-induced B str with zero 
average, and over fluctuations of electron density. These 
together determine the local values of the inter-LL sep- 
aration, ilo. a = \/2hc/lB a , and the local filling factors 
of the n — and n = ±1 LLs in each of the two val- 
leys a = K, K' . Eq.([T]) can also be obtained from Eq.Q 
neglecting strain and carrier density fluctuations. 

To measure the EPC, we fit Eq.© to the data of FigH 
The high-field anticrossing branches are most ac- 

curately resolved and therefore most suitable for fitting. 
The energy of the unperturbed electronic transition Qq is 
calculated using c = 1.08 x 10 6 m/s, which corresponds to 
the MPR at 25T. For low carrier density (0.4x 10 12 cm -2 ), 
we obtain g% + = 17cm" 1 /T 1 / 2 and gf = 12cm-7 Tl/2 - 
For intermediate carrier density (2.0 x 10 12 cm -2 ), fitting 
yields gfi + = 21cm- 1 /T 1/2 - Plugging the fitted g a± into 
Eq.([T]), we find v at the resonance B, and extract a car- 
rier density^ 0.4 x 10 12 and^ 1.6 x 10 12 cm~ 2 for low and 
intermediate density samples, respectively, in excellent 
agreement with those derived from the Raman spectra. 
Taken together, these results allow us to deduce the di- 
mensionless EPC. Notably, the Ar we obtain for the low 
and intermediate density samples differ by less than 5%. 
We get Ar ~ 0.028, in remarkable agreement with that 
measured from FWHM(G) in undoped samples at B=0, 
and predicted by density functional theory (a. M. [l3j. 

Finally, we discuss the observation of Raman scattering 
in the middle of the MPR anticrossing gap. The spectra 
for the intermediate carrier density show an additional 
component, indicated in blue in Fig(5pl. This is detected 
over a narrow B range, between 23 and 27T, reaching 
maximum intensity at 25T. Though we are not able to 
spectrally resolve this component at lower carrier densi- 
ties, a similar sharp increase of intensity is observed in 
the vicinity of the MPR at 25T (FigEJi, bottom panel). 
Eq. © shows that a finite Raman scattering inside the an- 
ticrossing gap is a specific signature of MPR in strained 
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FIG. 4. MPR in SLG with (p s ) = 1.6 x 10 12 cm 
presence of inhomogeneous strain, (-B 2 j r ) 1//2 = 6T 
lated using Eq.© with Sir ~ 1581cm" 1 , 7r « 12.3cm" 1 , 
7e = 50cm -1 , and (<5p 2 ) 1/2 = 3 x 10 n cm" 2 . Symbols are 
experimental results. The inset shows the calculated Ra- 
man intensity, neglecting LL broadening, depolarisation, in- 
homogeneity, but with strain-induced inhomogeneous B s t r - 
Dashed lines indicate the unperturbed E2 9 phonons and (0,1) 
magnet o- excit ons 



SLG. FigU] plots the Raman intensity calculated for in- 
termediate carrier density. To model the experiment, we 
assume density fluctuations described by a normal dis- 
tribution with a (p s ) ~ 1.6 x 10 12 cm -2 average and 20% 
standard deviation, imperfect light polarization leading 
to a 75%/25% mixture of a + and <7~ spectra, LL and 
phonon broadening of 7 e =50cm" 1 and 7r=12cm~ 1 
and, most importantly, (-B 2 ^) 1 / 2 = 6T. Away from the 
f2o = Sir MPR intersection, the Raman intensity follows 
the expected anticrossing MPR behavior, while the in- 
creased scattering in the middle of the anticrossing gap 
is due to the overall effect of inhomogeneous strain. This 
is further illustrated in the inset in Fig|4l with a Raman 
intensity map calculated for the ideal case of uniformly 
doped SLG (Sp s = 0) exposed to the same random B str , 
neglecting broadening (7r = 7e=0) and depolarization. 

In summary, we used polarization-resolved high-field 
magneto-Raman to investigate magncto-phonon reso- 
nances in graphenc. By varying the filling factor, we iden- 
tified different types of G-peak magnetic-field dependen- 
cies, providing a comprehensive experimental evidence of 
MPR on circularly polarized phonons. We also detected 



an unexpected increase of Raman intensity in the middle 
of the MPR anticrossing gap and assigned it to mixing 
and splitting of electron-phonon coupled modes caused 
by fluctuations of strain-induced pseudo-magnetic fields. 
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